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The dynamical structure factor of the frustrated spin 1/2 Heisenberg model for a series of frustration 
parameters at finite temperatures is presented. A sharp upper boundary of the spinon continuum 
is found. A simple method to extract the nearest neighbour coupling J from the spectral width is 
suggested. Signatures of frustration are discussed and a comparison with neutron inelastic scattering 
data is given. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.40.Mg, 75.50.Ee 
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The recent widespread interest in low-dimensional 
quantum spin systems has been brought up and main- 
tained by the discovery of new compounds with a variety 
of spin-spin couplings. Particular effort has been concen- 
trated on the first two inorganic substances undergoing 
a Spin-Peierls transition: CuGeC>3 and a'-NaV2 0s. 

There is extensive literature on the low temperature 
properties of CuGeC>3 (see for example Ref. |jj and ref- 
erences therein), but there is still no agreement on the 
model describing the uniform phase ( T > Tgp ) also. 
If the interchain interaction is negligible the most likely 
candidate is the isotropic frustrated spin 1/2 Heisenberg 
model |H in one dimension 



H = 2 J 



N 



Si ■ 5*2+1 + otSi ■ Oj+2 



(1) 



Alternatively Uhrig proposed a model including inter- 
chain interaction with a small value of inchain next to 
nearest neighbour coupling a Q] . 

In a recent study |^] it has been demonstrated com- 
pellingly that in the framework of the first model (Eq. 
(||)) for the unique set of parameters J/kp — 80K and 
a = 0.35 excellent agreement between experiment and 
theory is achieved for the magnetic susceptibility. 

In the following we focus on the dynamical properties 
of both substances at temperatures ksT > 0.6J. We 
compare our theoretical results with data obtained from 
neutron inelastic scattering (NIS) experiments Q per- 
formed on CuGeOa. There are still no corresponding re- 
sults for a'-NaA^Cv If it is confirmed that a'-NaV20s 
is well described by the pure isotropic Heisenberg model 
our results for a = can be taken as a prediction for this 
compound. 

We have determined the dynamic structure factor 
S(q, u>) at finite T for the frustrated Heisenberg model 



and the AAA-model for chains with N — 16, 18 
spins and periodic boundary conditions by the method 
of complete and exact diagonalization (CED) as de- 
scribed in detail in Ref. Q]. This was done for a — 
0.0,0.06,0.12,0.18,0.24,0.35,0.45,0.5. The largest sub- 
spaces to be treated are of dimension 415 for N = 16 and 
1367 for = 18 (see also Ref. ||). Our aim is to study 
the dependence of S(q, lj) on the frustration parameter 
a at T ^ and to highlight the unambiguous features 
distinguishing models with varying degree of frustration. 

We study chains with maximal 18 spins. Lanczos stud- 
ies of ground state properties of spin i systems routinely 
handle chains of 30 spins. So one could dismiss the 
present work as not worthwhile because it is expected 
to be too inaccurate. But thermodynamical observables 
can be determined by CED with high accuracy for tem- 
peratures ksT > 0.6J. An extensive investigation of 
finite size errors of dynamical correlations convinced us 
that this method, if carefully applied, gives sufficiently 
accurate results with easily controllable errors, which are 
not obscured by more or less arbitrary approximations. 
The obvious reason is that the number of contributing 
matrix elements is for T > many orders of magnitude 
larger than for T = 0. We refer to Ref. where a thor- 
ough finite size analysis is presented and to remarks in 
the course of this paper. With substantially more effort 
marginally longer chains could be studied exactly. But 
probably the additional accuracy would barely justify the 
effort. The real challenge is instead to extract as much 
reliable results as possible from the existing data. 

For ksT < 0.6J the method applied by us gives only 
poor results. We therefore restricted ourselves to temper- 
atures ksT > 0.6 J and considered qualitative and global 
features respectively. For higher temperatures where our 
method is accurate the existing experimental results are 
of qualitative nature only and our data can serve in the 
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FIG. 1. Dynamical structure factor for a = 0.0 (black 
solid line) and a = 0.35 (red dashed line), as a function of 
u/J for ksT = 3.8J. The triangles mark the upper boundary 
of the two-spinon continuum for the XXX-model. The blue 
vertical lines mark the experimental spectral boundaries read 
off from Fig.4 in Ref. §. 



future to disentangle magnetic and phonemic contribu- 
tions in the analysis of experimental data. 

We find that the structure function has for all practi- 
cal purposes (but not strictly) a bounded support in the 
q — uj plane. The lower boundary - well known from spin 
dynamics at T — - disappears gradually for growing T 
(see Figsjl], ||), whereas the upper boundary is nearly 
unchanged up to infinite temperature. The structure 
function develops for increasing T a steep slope at its 
boundary. We emphasize that the boundary as well as 
the functional form of S(q, u>) near the boundary deter- 
mined by CED show very small finite size errors: While 
the precise positions of local peaks of S(q, ui) for small 
and intermediate u are not perfectly iV-independent for 
small T, the shape near the boundary and the location 
of the boundary are remarkably stable. In any case the 
spectral boundary is easily determined independently of 
the criteria applied. This is in sharp contrast to the sit- 
uation at T = where the standard tool to study spin 
dynamics is the recursion method f|, which works well 
for low lying excitations but is much less accurate for 
higher ones. 

It is a remarkable and unexpected result that for T > 
the spectral density is bounded by a curve similar (but 



FIG. 2. As Fig.0 but for k B T = 0.6J. The triangles mark 
the upper and lower boundaries of the two-spinon continuum 
for the XXX-model. The blue vertical lines mark the experi- 
mental spectral boundaries read off from Fig.fb in Ref. 



not equal) to the two spinon boundary of S(q,u>) in the 
XXX-model detected by Midler et al at T = 0. For 
T > this was first observed in Ref. |7j . 

As for T = this results from a suppression of the ma- 
trix elements |(f 1 53(g) for large energy transfer. The 
density of allowed pairs of excitations however is not re- 
sponsible for the steep drop of S(q, u>) near the boundary 
since it is still sizable for values of the energy transfer far 
beyond the boundary of the two-spinon continuum. This 
is shown in Fig. || for a — 0.0, 0.35. 

The spinon is a rigorously defined notion only in the 
context of exactly solvable models (ESM) where Bethe- 
Ansatz techniques work [pT|-p7[. It is indispensable to 
understand the physics of spin dynamics at T = as first 
noted and discussed exhaustively in , see also |l8],|l9| . 
As we have just demonstrated, certain aspects of spin 
dynamics at T = in ESM and T > ( up to infinite 
temperature ) even for not exactly solvable models are 
strikingly similar. This makes it evident that the spinon 
is a much broader concept than originally anticipated. 

Our most conspicuous new result is shown in Figs. [I], 
|],|]: the spectral width is essentially independent of a 
for q = ir. This is of eminent importance for the analysis 
of experiments: it allows to determine the coupling J 
and its temperature dependence from neutron inelastic 
scattering without any further assumption on the model 



2 




0.0 2.0 4.0 6.0 8.0 10.0 12.0 14.0 16.0 18.0 
m/J 



FIG. 3. S(q = tt/2, w) for a = 0.0, 0.18, 0.24, 0.35, 0.45, 0.5 
at ksT = 3.2J and the density of contributing excitations for 
a — and a = 0.35 (rescaled). 




obtain the value of a from finer details of S(q, lu). These 
are: 

(i) The theoretical width of the spectral density for q ^ 
7r shows for low temperatures a distinct a-dependence. 
We extracted the boundary of the spectral density from 
the color contour maps of Ref. ||. They are shown as 
vertical lines in Fig. [I and ||. For T = 50K they coincide 
perfectly with the theoretical result valid for a = 0.35. 

(ii) The integrated intensity S oc (q) is shown in Fig.^ 
for T = 50K in conjunction with data from Arai et al. 
(see Fig. 2 of Ref. [||). If the vertical scale of the ex- 
perimental data is appropriately adjusted the agreement 
of experiment and theory is remarkable for a 0.35, 
whereas a = is clearly ruled out. There is no ques- 
tion of finite-size errors, they are much smaller than the 
experimental errors as shown in the figure. 

(iii) A further distinguishing qualitative feature is the 
trough-shaped structure detected experimentally which 
is fully developed at T = 300K (see Fig.4 of Ref. @). 
Fig. [| shows that the same phenomenon is present for 
a > 0.35 whereas for a = we find a hill-like structure. 

Summarizing we have found an unexpected high degree 
of correspondence of experimental data and theoretical 
results. We state that the model parameters J and a 
obtained from thermodynamical properties of CuGeC>3 
are also strongly favoured by NIS data. 

We conclude this section by predicting some properties 
of the NIS cross-section of a'-NaV20s. Assuming a = 
we find : 

1.) The spectral boundary wg±Awt,fl for fc^T = 1.2 J 
and q = mir/8 is given in the following table. Here Alor 
is much larger than Alol, because S(q,uj) for a = has 
a flat tail and becomes steeper for larger temperatures 
only. 



FIG. 4. S{q = tt,lo) for a 
and k B T = 3.2J. 



0.0,0.18,0.24,0.35,0.45,0.50 



parameters. 

A similar statement can be made about the normalized 
integrated intensity 



£(<?,w) = 



1 



Seo{q) Jo 



S(q, (jj')duj' 



(2) 



with S 00 {q) = Jq 00 S(q, ui')duj / . For q = n there exists 
for each temperature an interval were these functions are 
narrowly coalescing for all a values ( see e.g. Fig. ^|). 

Arai et al. J6| measured a spectral width of 32 meV at 
T = 5QK. This gives J/fcs = 79K and is in perfect agree- 
ment with J/ks = S0K determined from the observa- 
tion that the magnetic susceptibility has its maximum at 
T = 56K || . So we observe that completely independent 
experiments performed in the same energy range give the 
same value for the coupling J. At T = 300K the mea- 
sured width is 30 meV corresponding to J/ks — 

Having determined J by this very simple method, we 
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0.76 



2. ) The integrated spectral density Soo(q) at k R T = 
0.6 J is shown in Fig. g for a = 0. 

3. ) In clear distinction to CuGeC>3 the intensity is not 
suppressed but enhanced around ioj J = 3 and q = it. 

For a discussion of the underlying physics we refer to 
Ref. @. We add the following supplement. As we con- 
sider here solely isotropic models the spin S is conserved 
and we can study another classification of excitations 
contributing to S(q,w). There are two classes of tran- 
sitions behaving markedly different: 

1.) AS" = transitions: They do not occur for T = 0. 
In Fig. 7| it is shown that they contribute mainly to the 
central bulk of S(q,io). 
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FIG. 5. E(g = tv,u,T = oo) for N = 16 and a=0.0, 0.06, FIG. 7. S{q = tt,w) at fc s T = 3. 2 J for AT = 18 and 

0.12,0.18, 0.24, 0.35, 0.45, 0.5. a = 0.35. Contribution from transitions with AS = and 

AS = 1 
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FIG. 6. Soo(q) at k B T = 0.6J for a = 0.0,0.35,0.45,0.5 
for N = 16, 18 and experimental data at T = 50-K" read off 
from Fig. 2 in Ref. (The lower limit of integration for the 
theoretical curves is 3meV with J/ks = 80K in accordance 
with Ref. il.) 



2.) AS* = 1 transitions: Besides adding intensity to 
the central region they build up the steep slope near the 
spectral boundary. 

This note is exclusively devoted to an attempt to ex- 
plain NIS data of CuGe03 in the framework of model 
Eq. (|I|). We postpone other topics like a theoretically 
oriented study of spin dynamics in frustrated models as 
well as the application to e.g. Raman scattering to future 
publications. 
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